Theoretical de Haas-van Alphen Data and Plasma Frequencies of MgB 2 and TaB 2 
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The de Haas-van Alphen-frequencies as well as the effective masses for a magnetic field parallel to the crystallo- 
graphic c-axis are calculated within the local spin density approximation (LSD A) for MgB2 and TaB2. In addition, 
we analyze the plasma frequencies computed for each Fermi surface sheet. We find a large anisotropy of Fermi 
velocities in MgB2 in difference to the nearly isotropic behavior in TaB2. We compare calculations performed 
within the relativistic non-full potential augmented-spherical-wave (ASW) scheme and the scalar-relativistic full 
potential local orbital (FPLO) scheme. A significant dependence for small cross sections on the bandstructure 
method is found. The comparison with the first available experimental de Haas-van Alphen-data by Yelland et al. 
(Ref. 19) shows deviations from the electronic structure calculated within both L(S)DA approaches although the 
cross section predicted by FPLO are closer to the experimental data. The elucidation of the relevant many-body 
effects beyond the standard LDA is considered as a possible key problem to understand the superconductivity in 
MgB 2 . 
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I. INTRODUCTION 

The unexpected discovery of superconductivity in 
MgB2tJ, with a surprisingly high transition tempera- 
ture T c ~ 40 K, has raised considerable interest in 
the clarification of its electronic structure and its rela- 
tionship to the mechanism of superconductivity. The 
present theoretical approaches to this challenging prob- 
lem can be divided roughly into three scenarios em- 
phasizing: (i) a conventional electronic structure and a 
more or less standard electron-nhopon interaction within 
a (multiple) wide-bancLmodeln"!!!, (ii) the presence of 
strong electron-electronla ps magnetic intemctionsEl, and 
(iii) strong nonadiabatiola or polaronicO effects. In 
the present paper we shall mainly consider the first ap- 
proach which is based on calculations of the electronic 
structure within the local (spin) density approximation 
(L(S)DA) to the density-functional theory. The L(S)DA 
band structure calculationso revealed four sheets of the 
Fermi surface of MgB2 with a remarkable high anisotropy 
of Fermi velocities for c-holes. That anisotropy is impor- 
tant to explain the significant anisotropy of the upper 
critical field H C 2(T) in terms of a two-band Eliashberg 
modeH - LJ. The relatively high value of H C 2(0) is related 
to the strong electron-phonon interaction of the cr-holes 
on the tube-like Fermi surface sheetsQ. Details of the 
electronic structure are important for a precise under- 
standing of the pairing mechanism, like, e.g., by compar- 
ing the calculated "undressed" electronic effective mass 
with its actual dressed value due to the electron-phonon 
coupling. On the other hand, it is interesting to com- 
pare MgB 2 with related isostructural compounds having 
much lower transition temperatures or even showing no 
superconductivity at all. We choose here TaB 2 for which 
band structure calculations have been reported recently 



and for which single crystals can be producedO. 
It should be noted that stochiometric— ,TaB« is non- 
superconducting (at least down to 1.5 Kll3taO) and an- 
other phase, unidentified yet, should be responsible for 
the observed superconductivity with a T c of 9.5 Klia. 

II. METHODS 

Very detailed insight into the electronic structure 
can be gained from de Haas-van Alphen (dHvA) mea- 
surements. To the best of our knowledge there is as yet 
only one very recent measurement devoted to MgB2c2l 
and we are not aware of any dHvA-experiments for TaB2- 
With the aim to initiate possible experiments, we pro- 
vided calculated dHvA-frequencies and effective masses 
for the two materials under consideratior£3. From that 
data one can extract information about the Fermi veloc- 
ities averaged over the extremal orbits. We discuss also 
the possible enhancement of the effective masses by the 
electron-phonon coupling which is expected to be quite 
different for the two materials. Information about the 
Fermi velocities can furthermore be obtained by means 
of the plasma frequencies detectable from the Drude part 
of the optical conductivity. We compare also the plasma 
frequencies computed for each Fermi surface sheet and 
find a large amount of anisotropy in MgB2 in difference 
to TaB 2 . 

The band structure calculations have been per- 
formed by means of two computational codes: (i) 
the full*, relativistic augmented spherical waves method 
(ASW)IEl and (ii) the scalar relativistic full potential lo- 
cal orbital scheme (FPLO). In case (i) we used the von 
Barth-Hedin parametrization of the exchange-correlation 
potential. A fc-mesh of 148 special fc-points in the irre- 
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ducible Brillouin zone—part was chosen. For details of 
case (ii) see e.g. Ref.BE3. 

Both compounds crystallize in the hexagonal space 
group P6/mmm with lattice constants a — 3. 083 A, 
c = 3.5208 A and a = 3.082 A, c = 3.243 Alia for 
MgB2 and TaB2, respectively. The resulting ASW band 
structure J'or_MgB2 is similar with those published in the 
literature! ~ El and it is shown for the sake of complete- 
ness in Fig. 1. The relativistic band structure of TaB2 
is presented in Fig. 2. We have checked our results by 
comparing them with the scalar relativistic full-potential 
local orbitallFPLO) method repeating previous calcula- 
tions (MgB 2 ErQ, TaB 2 li3) and we found some differences 
for both compounds. In MgB2 these differences are not 
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FIG. 1. The calculated energy bands of MgE>2. The Fermi 
energy is at zero energy. 



due to the spin-orbit coupling in a fully relativistic cal- 
culation. Instead we are forced to ascribe these effects 
mainly to a slightly different crystal field potential in 
spherical or full-potential descriptions, which affects the 
relative position of boron derived n and a bands. The 
spin-orbit coupling of the Ta 5d states in TaB2 results in 
differences compared to scalar relativistic calculations, 
especially between T and A (compare Fig. 2 and Fig. 3 
of Ref. |l2|). Although these differences should not af- 
fect the absence of superconductivity in this compound, 
they lead to consequences for the Fermi surface to be dis- 
cussed below. To predict the dHvA-frequencies F for a 
magnetic field applied, say along the crystallographic c 
direction, we calculate the area A of any extremal orbit Q, 
according to the Onsager relation (written in cgs units) 
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It is convenient to measure F in units of the applied 
magnetic field and A in units of squared wave vectors, 
e.g., in kT (lkT=10 7 G) and the cross-section A in nm -2 , 
respectively. Then Eq. (1) reads simply 



F = 0.104867 x A 
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The effective mass of the corresponding orbit is given by 

h / dk de k 
m = — <b —r, with v k = -^r- 

27r/ |i; fe | dk 



(2) 



From the dHvA-frequency F and the effective mass m 
one can introduce an effective Fermi velocity 
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FIG. 2. The same as in Fig. 1 for TaB 2 



which corresponds to the hypothetical value for a cir- 
cular orbit with constant This expression is fre- 
quently used to estimate the Fermi velocity from dHvA 
data. Alternatively, one can also average 1/|?4| around 
the circumference and define 
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Obviously, both velocities defined by Eqs. (3) and (4) 
coincide for a circular orbit. For an orbit with an ideal 

0.9523 is estimated 



hexagon shape Vf/vf = y 7r/2\/3 
from Eqs. (2,3). Both of these quantities contain aver- 
ages of point properties of the Fermi surface around a 
cyclotron orbit. 

Information about the Fermi velocity is also con- 
tained in the plasma frequencies Lu p ^ a (a = x(y) or z) 



d 3 k 
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where we will distinguish the contributions from different 
Fermi surface sheets in the following (denoted by ujp a ). 
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Analogously to Eq. (5), components a of the Fermi ve- 
locities v a = y/ < >fs can be introduced which are 
averaged over selected sheets or over the whole Fermi 
surface (FS) 

/r] 3 k 
—S(e k -E f )(v k , a f. (6) 

The ratio of these Fermi velocity components related to 
the most anisotropic er-tubes with the strongest el-ph in- 
teraction at the same time can be used to describe upper 
bounds for the anisotropy parameter 7 of the upper crit- 
ical field H C 2, 



1 = KtlKl <\J< V 2 f>xy >tubes /2 < V 2 f z > tubes , 

if that anisotropy is the predominant one as suggested 
by several microscopic calculations^. 

III. MAGNESIUM DIBORIDE 

Let us first discuss the results for MgB2. The band 
structure near the Fermi energy Ef is depicted in Fig. 
1. The corresponding dHvA-data for a magnetic field 
directed along the oasis are shown in Table I. (Note that 
in our previous paperE3 the numbers for the orbits 2 and 
4 were mistakenly interchanged.) The Fermi surface of 




FIG. 3. (Color) The calculated Fermi surface, distribution 
of Fermi velocities, and extremal orbits (pink lines, the num- 
bers denote the orbits that are explained in the text) (mag- 
netic field || c-axis) of MgE>2. The Fermi velocities expressed 
by different colors (blue-slow, green-medium, read-fast) are 
given in units of 10 5 m/s as shown in the legend. 




FIG 4. (Color) The same as in Fig. 3 for TaB 2 . Note 
that for visibility purposes A is at the center of the hexagonal 
prism and the V point is the midpoint of the lower and the 
upper hexagon. 



MgB2 (see Fig. 3 ) consists of four sheets and there are 6 
extremal orbits for a magnetic field in c direction. We will 
denote them by numbers 1 to 3 centered around T and 4 
to 6 centered around A with increasing area in each case. 
So one can see that orbits 1 and 4 (2 and 5) belong to the 
smaller (larger) tube, where the remaining orbits 3 and 

6 are rather large and belong to the ring-like parts of the 
Fermi surface. The electrons on orbit 6 have the largest 
Fermi velocity and Vf of orbit 3 is smaller compared to 
that of orbit 6 due to the larger effective mass. The elec- 
trons on the tubes are slow in general, with the minimal 
velocity on orbit 2. The ratio of Fermi velocities for the 
two different groups of electrons is roughly 2. In addition 
there is a rather different strength of the electron-phonon 
interaction on both types of Fermi surface sheets: strong 
on the tubes and weak on 7r-derived rings. This confirms 
the introduction of a-jtwo-band model proposed for MgBj 
first by Shulga et alB and elaborated later on in Ref. p2. 

The relative difference (vf — Vf)/vf does not exceed 

7 percent and it is noteworthy only for those orbits which 
have a pronounced hexagonal form than a nearly circu- 
lar one (for example orbit 3 is such a aearly circular or- 
bit). According to the two-band modelQ we would expect 
a different electron-phonon mass renormalization (1+A) 
for electrons on the tubes (with A~ 1.2 •• ■ 1.5) and on the 
large orbits 3 and 6 (A~ 0.3). The comparison with the 
experimentally determined cross sections shows signifi- 
cant deviations from our L(S)DA predictions. The mi- 
croscopic reason for that unexpected deviation remains 
unclear at the moment. Standard electron-phonon in- 
teraction does not affect the cross-section. It might ex- 
plain only the renormalization of the Fermi velocities and 
of the masses. The cross-sections might be affected by 
strong electron-electron interaction and/or the presence 
of electronic instabilities competing with superconductiv- 
ity which might result in local changes (destruction) of 
the Fermi surface by a partial dielectrizatiqn, (gapping) 
from a charge density or spin density waveEia. The ex- 
periments yield a much narrower smaller tube than the 
L(S)DA calculation predict. While this already indicates 
a discrepancy between theory and experiment, more se- 
rious might be the fact that the wider, second tube is not 
observed in the experiment at all. Would this fact be con- 
firmed also for cleaner samples, i.e. the missing cross sec- 
tions with lower Fermi velocities could not be attributed 
to the residual impurity scatteringEj, it would imply the 
failure of the L(S)DA to describe correctly MgB 2 . 

We consider the elucidation of the responsible 
many-body effects which are missing in the standard 
LDA approach to MgB2, i.e. treating it as a simple sp- 
metal, as a possible key problem for our physical under- 
standing of the electronic structure of MgB2 and its re- 
markable superconductivity as well. Anyhow, in view of 
the great importance of the experimental dHvA data, the 
confirmation of a more or less unusual electronic struc- 
ture in MgB2 by other groups using different single crys- 
tals would be highly desirable. 

The experimental determination of the effective 
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masses by dHvA measurements and its comparison with 
the calculated electronic values gives important informa- 
tion on details of the electron-phonon coupling. If the 
assignment proposed in Ref. 19 is correct, strongly renor- 
malized Fermi velocities by a factor of 2 to 3 should be ex- 
plained. Comparing the calculated and measured masses 
for the smaller tube, one arrives at a slightly larger lo 
A w 1.27 compared with A ~ 1 obtained by Choi et al 
within full anisotropic (harmonic) Eliashberg theory. 

If the corresponding coupling constant for the larger 
tube would be also enhanced, possible contributions from 
the coupling of a soft mode should be considered. Oth- 
erwise it would be unclear, why for a non-enhanced 
Coulomb repulsion (measured by the Coulomb pseudopo- 
tential fj,*) or the absence of any other enhanced pair- 
breaking, the critical temperature is not higher than 40 
K. In this context the recent experimental observation of 
a soft-mode near 17 rpeV in Raman spectroscopy^, poini 
contact spectroscopy^ as well as in neutron scatterings 
is of considerable interest. Especially with respect to the 
assignment of that boson-mode, since standard harmonic 
phonon calculations are unable to explain the presence of 
such a low-lying modeu. Notably, also from recent stan- 
dard tunneling experimentsEH (unfortunately with a cut 
off at 26 meV, i.e. just slightly above the low-frequency 
anomalies) significant deviations from the standard har- 
monic phonon calculations have been reported. 

The anisotropy of Fermi velocities of MgB2 is even 
more visible in our calculated plasma frequencies (see Ta- 
ble II) , which are similar to the results of Reik. [|,| 28 2^ 
(see Table III). In our previous calculationE] a scaling 
factor of c/a was missing for U) px iy\. Here we have cor- 
rected the present results for that factor. It is another 
kind of anisotropy: the ^-component of Vf is about 6 
times smaller than the x(y) components for the two tubes 
(sheets 1 and 2). The anisotropy ratio v p tX t y )/ v p,z ( v a 
|w™ Q |) is much smaller for the sheets 3 and 4 and its 
average Fermi velocity is larger. The total plasma fre- 
quency, however, shows only a small anisotropy of 1.04 
between the z and x(y) components. In Table III we 
compare the anisotropy in the components of the total 
plasma frequency as obtained by various groups. The 
obtained numbers are generally in good agreement with 
one another. 

The anisotropy of Fermi velocities is much less pro- 
nounced in TaB2 (Tables III and IV) . Our fully relativis- 
tie_Fermi surface coincides almost with that published 
ir£3 and consists of three sheets. The extremal orbits 
are characterized as follows (s. Fig. 3): numbers 1 to 3 
are centered around the A point and correspond to the 
three different sheets. Orbit 4 is centered around K and 
orbit 5 corresponds to the small Fermi surface part in the 
neighborhood of K. There is one maximal orbit (number 
6) with nontrivial value k z (0.20615 x 2w/a). Since Ta 
is a 5c? element, the influence of the spin-orbit coupling 
is more pronounced in TaB2 than in MgBp. With re- 
spect to the scalar relativistic calculations^ we find two 



differences that affect the Fermi surface: First, a split- 
ting of bands along the high symmetry direction T — A, 
which leads to a lifting of band degeneracy at the Fermi 
level. Second, also at the if-point a small band shift 
occurs. The first spin-orbit modification of the Fermi 
surface makes one Fermi sheet disjointed along T — A. 
The calculated areas, effective masses and Fermi veloc- 
ities of TaB2 are collected in Table IV and one can see 
that any kind of anisotropies of Fermi velocities is much 
smaller than for MgB2. According to the calculation of 
the electron-phonon coupling constants^ we would ex- 
pect a rather small renormalization of effective masses in 
TaB 2 (A~ 0.05 • • • 0.2) again in difference to MgB 2 . The 
plasma frequencies are collected in Table V. Since there 
are no analogous FS sheets to the two tubes of MgB2 we 
do not obtain a corresponding anisotropy in Wp X t y \ with 
respect to a;™ z . This can be understood due to the loss 
of the quasi two-dimensionality of the B cr-states in TaB2 
due to the hybridization with the Ta 5d states. 

The hole Fermi sheet around A shown in the mid- 
dle panel of Fig. 4 resembles a feature calculated within 
an augmented-plane-wave code for ZrB2" The exper- 
imentally observed dHvA-frequencies of the central e- 
orbhxil (corresponding to orbit No. 2 in Fig. 4 | )- | near 
1.8 kT and 1.5 kT for the similar system TiB 2 E3, re- 
spectively, are significantly smaller than 3.48 kT in our 
FPLO-calculation for TaB2. The opposite behaviour is 
observed for the \i— orbit (orbit No. 6 in Fig. 4) where 
the experimental values for ZrB2 and TiB2 read 2.46 
kT and 3.44 kT compared with our FPLO calculated 
3.44 kT for TaB2. Further theoretical and experimental 
studies should be awaited to answer the important ques- 
tion whether there are significant deviations even from 
full-potential L(S)DA-calculations only for the peculiar 
MgB2 or also for other diborides with weak electron- 
phonon interaction. 



IV. SUMMARY AND CONCLUSIONS 

In summary, we have calculated the dHvA- 
frequencies and effective masses as well as the plasma 
frequencies of MgB2 and TaB2 which would be worth to 
be investigated experimentally in detail. First experi- 
mental de Haas- van Alphen data for MgB2 point to sig- 
nificant deviations. The comparison of the experimental 
masses with L(S)DA calculation could give important in- 
formation]- pa the electron-phonon coupling. According 
to Refs.EfElil3 we expect remarkable differences between 
MgB2 and TaB2: Large enhancement for one group of 
electrons in MgB2 and a small enhancement for the other 
electrons in MgB2 as well as in TaB2. We found a re- 
markable anisotropy of Fermi velocities in superconduct- 
ing MgB2 which is in difference to the more isotropic 
behavior of non-superconducting TaB2. 
Note added. After finishing the present second version of 
our work we have learnt on a closely related work by 



4 



Mazin and Kortusil who also stress the necessity of full- 
potential calculations for MgB2 and arrive essentially at 
very close numerical values for the cross sections and 
the masses under consideration (see Table I). Remain- 
ing slight numerical differences, comparison with other 
full-potential bandstructure codes, and somewhat differ- 
ent physical interpretation will be discussed clscwhcrca. 
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TABLE I. The LSDA de Haas-van Alphen data of MgB 2 (magnetic field || c— axis). Available experimental values according 
to the assignment of Ref. 19 have been added for the sake of comparison. Calculated quantitities without explicit label have 
been computed within the ASW scheme. 



orbit 


Fasw 


Ffplo 


Fexp 


A/Abz 


m/m 


TTXexp j TfXo 


Vf 


^f,exp 


Vf 




[kT] 


[kT] 


[KT] 


% 






[10 5 m/s] 


[10 b m/s] 


[10 b m/s] 


1 


0.99 


0.79 


0.54 


1.98 


0.25 


0.57 


7.9 


3.2 


8.5 


2 


2.16 


1.63 




4.3 


0.55 




5.4 




5.8 


3 


32.88 


34.2 




65.6 


1.70 




6.8 




7.2 


4 


2.31 


1.84 


1.53 


4.6 


0.31 


0.70 


9.8 


3.3 


10.5 


5 


4.31 


3.44 




8.6 


0.61 




6.8 




7.3 


6 


30.20 


30.6 




60.2 


0.92 




12. 




12.8 



G 



TABLE II. Squared plasma frequencies of MgB 2 in [eV 2 ]. 
One eV for the plasma frequencies corresponds to 80.4- km/s. 

FS sheet W L v x( y )/v z 



1 (tube) 


9.639 


0.242 


6.31 


2 (tube) 


9.444 


0.323 


5.40 


3 (ring) 


8.620 


22.053 


0.63 


4 (ring) 


23.042 


24.6698 


0.97 


All 


50.75 


47.290 


1.04 



TABLE III. Comparison of several calculated results for 
the two components of the total plasma frequency uj px and 
of MgB 2 (ineV). 

0Jp, x (y) to P ,z Reference 

7T2 6~88 this work 

7.21 6.87 Ref. 3 

7.02 6.68 Ref. 4 

7.13 6.72 Ref. 28 

7.04 6.77 Ref. 29 



TABLE IV. Calculated de Haas- van Alphen data of TaB 2 . 



orbit 


Fasw 


Ffplo 


Aasw 1 A-BZ 


m/m 


Vf 


Vf 




[KT] 


[KT] 






[km/s] 


[km/s] 


1 


2.08 


2.50 


0.041 


0.221 


1320 


1410 


2 


2.99 


3.48 


0.059 


0.356 


979 


1050 


3 


8.40 


8.31 


0.167 


0.677 


864 


926 


4 


13.26 


13.84 


0.264 


0.558 


1320 


1410 


5 


1.05 


1.00 


0.021 


0.236 


869 


891 


6 


1.15 


3.44 


0.023 


0.632 


339 


365 



TABLE V. Squared plasma frequencies ui px , ui pz , and the 
anisotropy of averaged Fermi velocities for TaB 2 (in eV 2 ). 



FS sheet 






vx(y)/v z 


1 


68.125 


84.779 


0.90 


2 


15.149 


17.049 


0.94 


3 


3.243 


3.254 


1.00 


All 


86.52 


105.08 


0.91 
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